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Abstract 

A method is presented in which matrix elements for some processes are calculated 
recursively. This recursive calculational technique is based on the method of basis spinors. 
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1 Introduction 

o 

The possibility of investigation higher and higher energies at present and future colliders, entails 
the necessity of predicting and calculating with high precision more and more complicated 
processes. When the number of final particles is high it becomes hard, even to calculate the 
corresponding tree level Feynman diagrams and the final expression of cross section is often an 
intricated function of several variables, inadequate for practical use. This has lead to that it is 
necessary abandon the standard methods for perturbative calculations and to use instead the 
\jT) ', new effective ones. 

The standard method to obtain a cross section with the fermions in perturbative quantum 
' field theories is to reduce the squared amplitude to a trace from products of 7-matrices. An 

alternative approach is to calculate the Feynman amplitudes directly. The idea of calculating 
amplitudes has a long enough history. In 1949 it was suggested in Ref. P to calculate a matrix 
element by means of explicit form of 7-matrices and Dirac spinors (more detailed bibliography 
Ph| on the problem can be found in Refs. [2ED- 

Various methods of calculating the reaction amplitudes with fermions have been proposed 
and successfully applied in recent years. In general the methods of matrix element calculation 
can be classified into two basic types. The first type includes methods of direct numerical cal- 
culation of the Feynman diagrams. The second type includes methods of analytical calculations 
of amplitudes with the subsequent numerical calculations of cross sections. Notice that there 
are methods of calculating cross sections without the Feynman diagrams jlJ|S]. 

Analytical methods of calculating the Feynman amplitudes can be divided into two basic 
groups. The first group involves the analytical methods that reduce the calculation of S'-matrix 
element to a trace calculation. The reduction of a matrix element to trace calculation from 
products of 7-matrices underlies a large number of methods (see, Refs. [21 El and Refs. [U]-[T3] 
etc.). In this method the matrix element is expressed as the algebraic function in terms of 
scalar products of four-vectors and their contractions with the Levi-Civita tensor. 

The second group involves the analytical methods that practically do not use the operations 
with traces from products of 7-matrices. The method of the CALCUL group which was used 
for the calculations of the reactions with massless fermions is the most famous among [Tij-|18j. 
The basic idea behind the CALCUL method is to replace S'-matrix element by spinor products 
of bispinors and to use the fact that expressions u\ (p) U-\ (k) are simple scalar functions of 
the momenta p, k and the helicity A. However, the operation of matrix element reduction is 
not so simple as the calculation of traces. It requires the use of Chisholm spinor identities 
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(see [TED- Also it takes the representation of contraction fa = p^^i^ with four- momenta p^ 
and polarization vectors of external photons through bispinors. For gauge massive bosons the 
additional mathematical constructions are needed [16J. 

There are generalizations of the CALCUL method for massive Dirac particles both for special 
choices of the fermion polarization (^H],PH|-[22|) and for an arbitrary fermion polarization 
121] • We call the polarization states of fermions in Ref-O EE] as Berends-Daverveldt- 
Kleiss-Stirling or BDKS-states. 

Notice that for massless fermions we can obtain amplitude in terms of the scalar products 
of four-momentum vectors and current-like constructions of the type J 7 ' ~ U\ (p) 7 m ^a {k). The 
components of J M are calculated by means of momentum components p, k (so called i?-vector 
formalism; see Ref. [2*o]). 

For BDK S-states Ref. [26J presents the iterative scheme of calculation that reduces expres- 
sion for the fermion chain u\ (p) Q u\ (k) to the combination of spinor products u\ (p) u\ (k) and 
(or) u\ (p) 7 M (gy + fiUTs) (k) by means of inserting the complete set of non-physical states 
of bispinors (with p 2 < 0) into the fermion chain. 

In all above-mentioned methods the spinor products and current constructions were calcu- 
lated by means of traces and then used as scalar functions of the momenta and of the helicities 
(similar to scalar product of four- vectors) . 

Due to their easy implementation methods of matrix element calculation have become a 
basis for modern programs dealing with evaluation cross-sections for various processes. Ex- 
amples for such programs are generators HELAS [27|, GRACE j2H], MadGraph ' MEGA [50] . 
FeynArts/FormCalc4 H7TFT21 and CompHEP [33J (planned the calculation of the matrix element 
|34|). There is large number of more specialized programs as AMEGIC++ (35], ALPGEN [36J, WPHACT 
|3*7j . LUSIFER jHH] and et al. The detailed list of such programs can be found in Ref. [3U] . 

In the paper we describe an approach to Feynman diagrams which is based on the using 
of an isotropic tetrad in Minkowski space and basis spinors connected with it (see |4(Jj.|41j). 
Here we don't use an explicit form of Dirac spinors and 7-matrices and the operation of trace 
calculations. The method is based on the active using of the massless basis spinors connected 
with isotropic tetrad vectors and we will call it as Method of Basis Spinors (MBS) . 

In this method as well as in the trace methods the matrix element of Feynman amplitudes 
is reduced to the combination of scalar products of momenta and polarization vectors. Unlike 
spinor technique in different variants Jlj-jStX, this method doesn't use either Chisholm identi- 
ties, or the presentation of the contraction p 1 with four vector p and of the polarization vector 
of bosons through the bispinors. Unlike WvD spinor technique JH],j2I] MBS doesn't use special 
Feynman rules for calculating of the matrix elements. 

We propose to use recursion relations as a technique to evaluate the Feynman amplitudes 
of processes. The advantage of recursive technique is that for calculation of a n + 1 matrix 
element of some process we can use the calculation of n process. Both for analytic and numerical 
evaluation this is asset. 
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2 Method of Basis Spinors 



When evaluating a Feynman amplitude involving the fermions, the amplitude is expressed as 
the sum of terms that have the form 

Mx p ,\ k (p,s p , k, s k ;Q) = 

= M Xp ,x k ([p] , [k] ; Q) = u Xp (p, s p ) Q u Xk (k, s k ) , (1) 

where X p and X k are the polarizations of the external particles with four-momentum p, k and 
arbitrary polarization vectors s p , s^. The operator Q is the sum of products of Dirac 7-matrices. 
The matrix element © with Dirac spinors is a scalar function. Thus, it should be expressible 
in terms of scalar functions formed from the spin and momentum four-vectors of the fermions, 
including p, s p , k, s k and the operator Q. 

We will now consider that in the our approach this matrix element © can be represented 
as linear combinations of the products of the lower-order matrix elements. 



2.1 Isotropic tetrad 

We use the metric and matrix convention as in the book by Bjorken and Drell j32], i-e. the 
Levi-Civita tensor is determined as e i23 = 1 and the matrix 75 = i7°7 1 7 2 7 3 . Let us introduce 
the orthonormal four-vector basis in Minkowski space which satisfies the relations: 

F ■ l» -V> 1% ■ 11 - q ■ % = <T , (U ■ Ib) = 9AB, (2) 

where g is the Lorentz metric tensor. 

With the help of the basis vectors I a (A = 0, 1, 2, 3) we can define lightlike vectors, which 
form the isotropic tetrad in Minkowski space (see, 



b P = — 2 — ,nx = 2 >(/ ? ' A = ±1 )- ( 3 ) 



From Eqs. ©, © it follows that 

(VM = %, (n A • n-p) = ^ , (b P -n x ) = 0, ( I ) 



5\, p , x _ 5\, P 



9^ = 2J2 [K ■ b- x + K ■ n-J . (5) 

A=-l 

It is always possible to construct the basis of an isotropic tetrad © as numerical four-vectors 

(6 ±1 ) M = (1/2) {1, 0, 0, ±1} , (n±0„ = (1/2) {0, ±1, i, 0} (6) 

or by means of the physical vectors for reaction. 

For practical applications it is convenient to introduce additional four-vectors 

b p = 2 bp, h x = 2n x (7) 



3 



By means of the isotropic tetrad vectors we can determine the polarization vectors of mass- 
less (also and massive) vector bosons. For photons with momentum k^ and helicity A = ±1 we 
use the following definition of polarizations in the axial gauge 

provided that, the four- vectors are linearly independent. 

2.2 Massless basis spinors 

By means of the isotropic tetrad vectors (JHJ) we define massless basis spinors u\ and u\ (pi) 

JLiua = , m a(M =^i«_a(6-i) , (9) 

u\u x (6±i) = u A (6±i) (10) 
with the matrix = 1/2 (1 + A75) and the normalization condition 

ux(b±i) u\(b±i) = u)\ ]l ±1 . (11) 

The relative phase between basis spinors with different helicity is given by 

flxu-u (b-i) = 6x,uU\ . (12) 

The important property of basis spinors Q is the completeness relation 

1 

£ wa(6 a )«-a(6-a) = /, (13) 

A,A=-1 

which follows from Eqs. (jSJ , p2jl . Thus, the arbitrary Dirac spinor can be decomposed in terms 
of basis spinors u\ (6a)- 

2.3 Dirac spinors and basis spinors 

Arbitrary Dirac spinor can be determined through the basis spinor @ with the help of projec- 
tion operators u Xp (p, s p ) u Xp (p, s p ). 

The Dirac spinors wf (p, s p ) for massive fermion and antifermion with four-momentum 
p (p 2 = nip) , arbitrary polarization vector s p and spin number A = ±1 can be obtained 
with the help of basis spinors by means of equation: 

At \ (P' + Arrip) (1 + A75 Q 
2\/( & -i " {P + rripSp)) 



Tx (p, S p ) U_AxA (6-1) • (14) 
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The notation (p, s p ) stands for either u\ p (p, s p ) (bispinor of fermion; A = +1) or v\ p (p, s p ) 
(bispinor of antifermion; A = —1). Here we have introduced the abbreviations 

& = P - ±J Y^ ■ (15) 

The bispinors u\ (p, s p ) and v\ (p, s p ) satisfy Dirac equations and spin conditions for massive 
fermion and antifermion 

i> U\ (p, s p ) = m p U\ (p, s p ) , i> v x (p, s p ) = -nip v x (p, s p ) , 
75 ^ P u x (p, Sp) = A M A (p, s p ) , 75 fS p v x (p, s p ) = A v x (p, s p ) . (16) 

We also found, that the Dirac spinors of fermions and antifermions are related by 

V\ (P, S p ) = -A75 U—X (p, s p ) , (P, Sp) = u_ x (p, s p ) A 75 . (17) 

Let us consider the particular case of Eq. (jl4|) -the BDKS polarization states of fermions, as 
they are the most-used in calculations of matrix elements. The polarization vector of BDKS- 
states is defined as follows [1311111211221: 

V mpb-x 

sks = Sp- ■ ' 1 ^ 1 



nip (p • ' 

Performing a calculation for Eqs.(|51). ()14Jl and Eq. (jl8|) . we find the simple result for massive 
Dirac spinor [1211231211: 

w x (p, s KS ) = 1 (19) 

(p ■ 6_i 



Notice that, in Ref.|26j relation between the Dirac spinor of fermion and the Dirac spinor of 
antifermion differs from the Eq. (jl7|) . 

The Dirac spinor ux (p) of massless fermion with momentum p (p 2 = 0, (p • 7^ 0) and 
helicity A is defined by (for example, see Ref.[T6j) 

u M = ^M. (20) 



(p ■ 6_] 



2.4 Main equations of MBS 

The spinor products of massless basis spinors (0) are determined by 

ux (b c ) u p (b A ) = 5 X ,- P S c ,-a, C,A = ±1, A,p = ±l. (21) 

With the help of Eq. (0) Dirac matrix 7^ can be rewritten as 

1 



1" 



[V-x%+ ■ (22) 



A=-l 
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Using the Eqs.(|9|).([Hl|) and (|T2"j) we can obtain that 

-fu x (b A ) = V> A U—\ (6_ A ) - A n^ AxA M _ A (6 A ) (23) 

and 

75 m p (6a) = pu p (b A ) . (24) 

Ens. (|23|l - (|2^l and Eq.flST]) underlie the method of basis spinors (MBS) |4U[ 14 1|. 

By means of the Eq. (|2H|) we can determine that product of two 7- matrices can be represented 

as 

YYux (6a) = Yjtf u x (b A ) - A X% x u x (6_ A ) , (25) 
where X^", are the Lorentz tensors: 

X A, A = K ■ K-Ax\ - ™-Ax A ■ K , ( 26 ) 

YZ u x =bt A -b A + ^ AxX .^__ AxX . (27) 
The product S n = 7^ 7^ . . . 7^™ can be written as 
S n u x (b A ) 

= Btx'^ u K (b K ) - A N { /X^ } u K (b^J , (28) 

where 

\' n = (-l) n \ , A' n = (-l) n A (29) 

and B^ 1 ^'"^, j^fjvi,-^} arg gome L oren tz tensors, which are related with isotropic tetrad 
vectors (J7J). 

As follows from Eqs.(J2HJ),(J23J) we have that in particular cases: 

4 W a } = 6>\ J\/ft>=n^ xA , (30) 

A, A A, A J A, A A, A ■ 

3 Recursion relations for the matrix elements 

The basic idea of Method of Basis Spinors is to replace Dirac spinors in Eq.(JJ) by massless 
basis spinors ©, and to use the Eq. (}2T]) and Eqs. p3*j) - (j2^j) for calculation of matrix element 
(0) in terms of scalar functions B,J\f. With the help of the Eq.(|14j) the matrix element (JTJ 
transforms to fermion "string" with massless basis spinors u x (b A ) i.e. 

M Xp ,\ k (p, s p , k,s k ;Q) = 

= u_ Xp (6_i) T Xp (p, s p ) Q T Xk (k, s k ) u_ Xk (6_i) = 

= -M_ V - Afc (6-1, 6_i ; T Ap (p, Sp ) Q T Xk (k, s k )) , (32) 

where operator T A is determined by Eq.(|14j). Let us consider a special variants of a matrix 
element. 
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3.1 Basic matrix element 

Let us consider an important type of a matrix element when p = b^c and k = h^, i-e. 

M-„, p Q>-c b A ;Q)=T^ p A [Q} 

= U— a (fe_c) Q Up (b A ) . (33) 

We call this type of matrix element as basic matrix element. 

Note that, the matrix element (fljl is a particular case of basic matrix element i.e. 

M Xp ,\ k (p, s p , k, s k ;Q) = T);~\ k [T Xp (p, s p ) Q T Xk (k, s k )] . (34) 

With the help of the completeness relation (J13|) we can obtain the recursion formula for 

r£ P A [QM 

r% A [QM= £ if^ [Qi] r°; p A [Q 2 ] . (35) 

-D,A=-1 

By means of the relations (|25 j) .(|25 j) .(p5 )l and Eq. (J2T|) it is easy to calculate in terms of 
the isotropic tetrad vectors. For instance, 

r? P A W = 5„, p (do, -aK-A 5 ca nt Axp ) , (36) 

r? f [V 71 = <W (5 c ,aY£:; - A 5c,-aX!£ v p ) (37) 

and 

= S 9t P , n (6 C , A' n B { A y-^ - A 8 C ,-A> n N { /X^) . (38) 

With the help of the Eqs. (J35|) and (|38|) we obtain recursion relations for g^'"' and 



A, A —A' , , A' t A, A + I i J ■ /V A' ,, A' t A, A ' 



The recursion Eqs. (|39j) - (j4(Jj) allow to convert scalar functions B,Af into Lorentz tensors 
in terms of isotropic tetrad vector with the help of Eq.(|3U|) or Eq.(|31|). For example, the 

constructions q^q^q^B^ 2 ' ^ = B { A q ^ q2 ' q3} and q£ 1 q% 2 qg a M]&' ^ ^ = ^aX^'^ can be 
represented as a combination of the scalar functions X, Y and scalar products of the isotropic 
tetrad vectors: 

9l H2 % °A, A — D A, A 

93 • &a) + (?3 ■ n-AxA) X^ 2 A , (41) 

„Ml „A»2_A»3 A/"{Mi> M2, M3} _ 92, 93> 

«1 «2 ?3 ■ /V A, A — - /V A, A 

= (q 3 ■ n-Axx) YX)> q l - (q 3 • &a) ^-a . (42) 
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3.2 Decomposition coefficients 

The next type of lower-order matrix element is 



M p , Xp (b A ,{p};I)=M p , Xp (b A ,[p]) 
= u p (b A ) u Xp (p, s p ) . 



(43) 



The matrix element ()43|) is determined by the decomposition coefficients of an arbitrary Dirac 
spinor u Xp (p, s p ) on basis spinors Q- 

With the help of the Eqs. ()14Jl . (J2UJl the matrix element ()43|) transforms to 



M p , Xp {b A ,\p}) = 

_ up (b A ) [f 1+ f_ x u_ Xp (&_, 

b-i ■ £? 

for the massive fermions with arbitrary vector of polarization and transforms to 

Up(b A ) i>u_ Xp (6_i) 



(44) 



M p , Xp (b A ,p) 



(45) 



b-x -p 



for massless fermions. 

Using the Eqs. (J23j) - (J2H) and ()25jl the matrix element (j4*3^) is reduced to an algebraic expres- 
sion in terms of scalar products of isotropic tetrad vectors and physical vectors or in terms of 
components of four-vectors. 

Let us consider massless fermions. Using Eqs. (J21|) .(j23 |) we obtain, that 



/ 



M P)X {b A ,p) = 6. 



A-p 



p ■ b-i ) + 5 A 



\ 



(p ■ n_ x ) 



p ■ b. 



\ 



'J J 



(46) 



For numerical calculations, as well as in the case of spinor techniques, it is convenient to 
determine the (J46j) through the momentum components p = (p° , p x = p° sin 9 p sin (p p , p y = 
p° sin 6 p cos ip p , p z = p° cos P ) 



M p> a (b A , p) = 5 X _ p 5 A -i\/p~ - S A>1 Aexp (-i\<p p ) a/p+ 



5 A -\ sin y - 5 AA A cos ^ exp {-iX<p p ) 



(47) 



where 



p ± =p°± p z , p x + i\p y = \J {p x ) + (py) exp (i\(f p ) • 



Let us consider massive Dirac particles with arbitrary polarization vector s p . After evalua- 
tions we obtain, that the decomposition coefficients for a massive fermion with momentum p, 
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an arbitrary polarization vector s p and mass m p can be written as scalar products of tetrad 
and physical vectors 



1 



6-1 ■ 



*x, P K-i (6-1 " en + 5a,i (n-x p ■ £?) } + 



M p , Xp (b A ,p, s p 



+o\, p j<u,i y -i, - Xp + ^^ A -i,-Ap 

where the scalar functions Y p ' q ,X p ' q are determined by Eqs.(|2fi )) - (j27j) . 

For BDKS polarization states with the polarization vector (JT5J) the matrix element 
has a compact form 



(48) 



M Pi \(b A , p, s KS ) = S x , —f_ 



$a,-i\ [p- 6_i ) + 5 A , 



(p ■ n_ 



p ■ b\ 



+ 



+5\, P 8a,i 



The matrix element 
Eq.(dZD: 



p ■ b\ 



(49) 



with the antifermion can be easily obtained with the help of 



M P} Xp (b A , [p]) = u p (b A ) v Xp {p 
= p\ p M Pt -x p (b A 



(50) 



3.3 Recursion relation 

With the help of completeness relation (|T3|) the amplitude (JTJ) with Q = Q2Q1 is expressed as 
combinations of the lower-order matrix element 



Mx p ,\ k 



[k};Q 2 Qi) = 
M Xp , a ([p] , b A ; Q 2 ) M- a , x k (&-a, [k] ; Qi 



(51) 



(T,A=-1 



This insertion allows us to "cut" fermion chain into pieces of fermion chains with basis spinors 
u\(b A ). Hence our formalism enables to calculate the blocks of the Feynman diagrams and 
then to use them in the calculation as scalar functions. All possible Feynman amplitudes can 
be built up from a set of "building" blocks. 

Let us consider the matrix element (0) with an operator 



Z^ 1 ' — Q n Q n -i ■ ■ ■ Q\Qo 



(52) 
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with Q = I. In the Eq.(|52|) all operators Qj have an identical mathematical expressions. Using 
Eq. flgID we find that 

M Xp ,, k ([p],[k];Z^)^M^ Xk ([ P ],[k]) = 
i 

= E M^ t0 (\p)M)M™^{b- A ,[k\) , (53) 

a,A=-l 

where matrix element .M^ Afc [&]) can be calculated with the help of recursion relation 

a<SL ( 6 -a,[*]) = 

= E r £f M -p,1 ( b ~c> M) ■ (54) 

p,C=-l 

Once scalar functions .M Pi \ k (b A , [k]) and [Qj] are known (see EqB.(|3Sj). and 

Eqs. (|S^l - (|S7jl ). it is possible to evaluate the higher order of M._ aX [k]) with the help of 

recursion relation (J54j) . 

4 Examples 

Consider the "toy" example 

M"i \ k (p. fc > s fc) = 

= u Xp (p, s p ) i n . . . { x u Xk (k, s k ) , (55) 

where qj are some arbitrary four-vectors. 
Therefore, we have that 

Z^=i n i n _ x ...i x . (56) 
Using the Eqs.(j2SD-(EZD, (El) we find that 



M% Xk (6_c, [*]) = (fc ■ 6_c) -M^ ^ - 



p, *k 

and 



-C {qj-nojMWib-cM (57) 



a, (6-c, [*]) = ficj^ fe, . . . qj - k ) M^l Xk (b^_ h , [k]) + 
where CJ_ fc = (-l)*"* C, p$_ fc = (-l)'"* p, fc < j. 
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With the help of Eq. (|53p we can obtain the recursion formulas for calculation matrix element 



1 



= E «v 

p,C=-l 

-G ((&-n 0p )A<^ ) (6- ,[*])' 



(59) 



and 



1 



P,C=-1 



[*] + 



(60) 



We have obtained that the matrix element (jB3|) can be represented as a combination of the scalar 
functions B,Af, decomposition coefficients M.\ P , P (\p) , be) and lower order matrix elements. 

4.1 Photon emission 

Let us consider the matrix element where photon with momentum k and helicity a — ±1 
emission occurs from the incoming electron 



M X2 , Ai (M , M ; Q z k ) = u\ 2 (p 2 , s P2 ) Q z fc m Ai (pj, s pi ] 

Z k = e 



with 



(pi — A;) 2 — m 2 

Using the algebra of 7-matrix and Dirac equation the operator Z/, can be rewritten as 

(pi e a (&)) ft 4<j ( k ) 



Zu = —e 



Now we get 



( Pl k) 2 (pi fc) 
(Pi £ ff (A;)) 



(61) 
(62) 

(63) 



(pi fc) 



2(pi fc)' 

The recursion technique (see Eq. (l5T|l ) imply 



(64) 



A^a 2 , cr ([p 2 ] , Q) (b-c, N ; # ^ (k)) = 

C=-l 

1 



(65) 



-,c=-i 
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where 

(p, 8p) = 1 (66) 

(*-•«) 

for fermion with arbitrary polarization vector (see Eq.(|14|)). 

ft (p, s„) = ^±^- (67) 

(p • Li) 

for fermion with BDKS polarization vector (j!9j) . 

r A (p, Sl ,) = — =lL= (68) 

for massless fermion (see Eq. ()20p ). 

Let us consider initial massive fermion with BDKS polarization vector in expression (|65p. 
After calculations, for the matrix element 1)61)) with the BDKS polarization state of initial 
fermion and helicity a of photon we have the exact formula in terms of lower-order matrix 
elements with operator Q and scalar functions B, M: 

M X2 , Xl flpa] , [px] ; Q Z fc ) = -e ^M-^, a, (N , M 5 Q) + 

(pi //) 



+ 



A<A 8 ,AxW,6-i;Q)^-i;-rf* ) ' ,tt + 



2 (Pi fe) J (pi 6_i) 



+^a 2) a, ([p 2 ] , 61; g) P1 + m pi { m X2 , _ Al ([p 2 ] , 6_ i; g) K.t_X w 
+A^A 2) -A 1 (N,&i;g)x_ fc 1 ';f ) }] . (69) 

Scalar functions g fe ' ^W' 3 ^ jy^ fc > £ <>-w»pi are determined in terms of scalar product with the help 

of the Eqs.@]}-@2J and scalar functions Xj^j* , K_J'_% (fc) are determined by Eqs. (|25jl - (|27jl . 

Scalar functions can be easily calculated in terms of physical vector components. Using a 
vector of polarization of a photon as Eq.fjHJ) we obtain simple result: 

B\ L'W' p = V2 a {p-k T , a, [K, + k~p T , , 

Y_l ll[ k) = -y/2 Ax <f Alj _ CT A; T , _ CT , if ] = -y/2 5 Xl , -a k~ , (70) 

where 

p ± =p°±p z , p T , A =p x + i Ap y . (71) 
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4.2 The process e + e — > 

Consider the process 

e + (p 2 , cr 2 ) + e" (p u cri) -> 7 (fei, Ai) + 7 (fc 2 , A 2 ) H h 7 A„) , 

where the momenta of the particles and spin numbers are given between parentheses. 
The Feynman diagrams of the processes (J72)) contain the matrix element 

M (x^x 2 ,...x n) ^ Sp2 ^ i5 ^. ^ kn) = 
= MW, ([p 2 ] , ^]) = v a2 (p 2 , s P2 ) ? Xn (K) ■ ■ ■ 

• ' ' ?As (Ks) n2 _ o A 2 («2) n2 _^ 2 Ai «<n (Pi, Spi) + 

+ (n! — 1) other permutations of (1,2, ... ,n) , 

where 

3 

1=1 

Hence, we have that 

ryin) ®n + m Ql+m 



and 



Here 
and 



1 

= J2 •M CT2 ,p([P2],6c;A„(^))A^_ P)CT1 (6_ c ,N;^ (n - 1) ) - 

P,C=-1 



s p 2 ) 4x n (k n ) tip (be) 
Mt'H [Pi]) = {b-c, [Pi] ; ^ (n_1) ) • 



Using the expressions ()37)) and (JHJ) we obtain, that ()77)1 for BDKS massive Dirac 
(fTTJj) is determined by 

•Hr 2) p M > &C! (fen)) = 



= (». • 6-1) [<w + 

,1 (^cx 2 ' ^n) I 



with e n = e Xn (K 
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The final recursion relation of the process ()72|) with arbitrary helicities of the photons and 
BDKS polarization states of positron is written as 



M% ai (N , N) = vfip^i 

X P - 2 CM-^l Ox, N) + Yi P TMt-% (6-1, N) 

e n ) M%£ (6-1, tPi]) + (6-1 • en) (6 1; [ Pl ])) 



+m (n 



iO) 



where (see Eg. ([59)1 and Eq.tJHDJ)) the matrix element A^l^. (6_c, \pi]) with arbitrary polariza- 
tion vector of electron is calculated by means of the recursion formula 



(6-c, N) 



1 



Q] - rrv 



{ 



m [ (e j ■ b. c ) M^l\ (6_ c , [pi]) - C ( £j ■ n Cp ) A^J (6c, [*])] + 
+Ycf 3 M^ {bci M + G ^VM^ (6_c, tPi]) } • 



5 Summary and Acknowledgements 

In present paper we have formulated a new effective method to calculate the Feynman am- 
plitudes for various processes with fermions of arbitrary polarizations. In our method it is 
much easier to keep track of partial results and to set up recursive schemes of evaluation which 
compute and store for later use sub diagrams of increasing size and complexity. 
In our approach of the matrix element calculation: 

1 we don't use an explicit form of Dirac spinors and 7-matrices 

2 we don't use the calculation of traces 

3 as well as in the trace methods the matrix element of Feynman amplitudes is reduced to the 

combination of scalar products of momenta and polarization vectors. 

4 Unlike spinor technique in different variants in this method we don't use either 

Chisholm identities, or the presentation of the contraction p 1 with four vector p and of 
the polarization vector of bosons through the Dirac spinors. 

5 Unlike WvD technique JE],[2I] m this method we don't use special Feynman rules for calcu- 

lating of the matrix elements. 

6 Expression for matrix element calculates M.\ \ k (p,s p , k,Sk ;Q) (0) f° r a U values \ p , Xk si- 

multaneously. 

The recursive algorithms can be easily realized in the various systems of symbolic calculation 
(Mathematica, Maple, Reduce, Form) and in such packages as FeynArts [3T], FeynCalc ITf]. 
HIP and so on. 

I would like to thank organizers for warm an kind hospitality throughout the Conference. 
Also I want to thank A.L. Bondarev for his useful remarks and discussion. 
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